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Abstract 

Using in-situ synchrotron tomography, we investigate the coarsening dynamics of barium borosilicate 
melts during phase separation. The 3-D geometry of the two interconnected phases is determined thanks to 
image processing. We observe a linear growth of the size of domains with time, at odds with the sublinear 
diffusive growth usually observed in phase-separating glasses or alloys. Such linear coarsening is attributed 
to viscous flow inside the bicontinuous phases, and quantitative measurements show that the growth rate 
is well explained by the ratio of surface tension over viscosity. The geometry of the domains is shown 
to be statistically similar at different times, provided that the microstructure is rescaled by the average 
domain size. Complementary experiments on melts with a droplet morphology demonstrate that viscous 
flow prevails over diffusion in the large range of domain sizes measured in our experiments (1-80 ^m). 
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1. Introduction 

Liquid-liquid phase separation is an appealing 
mechanism for building materials with tailored mi¬ 
crostructure [T]. Molten materials such as poly¬ 
mers P], silicate melts PI a [5] or metallic al¬ 
loys iniiziiH] exhibit liquid-liquid immiscibility un¬ 
der a critical temperature [S]. During a quench 
under the immiscibility dome, spinodal decomposi¬ 
tion or nucleation and growth are first responsible 
for the growth of composition fluctuations. Then 


^corresponding author 


the size of domains of fixed composition increases 
in order to minimize interfacial energy - a later 
stage called coarsening HO]. The ability to tune 
the typical size of the microstructure by controlling 
coarsening kinetics is paramount for applications of 
phase separation, among which porous glasses m 
and membranes [12] , cellular materials for the food 
industry [Hill], core-shell particles quantum 
dots m or super-hydrophobic thin films m- 

For most materials, coarsening results from 
molecular diffusion through interfaces between the 
phases, transporting matter from highly-curved in¬ 
terfaces to flatter interfaces. This mechanism oper- 
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ates whatever the topology of the phases - either for 
interconnected domains or for droplets [18]. Dif¬ 
fusive coarsening is a slow transport mechanism, 
with the typical domain size growing with time t 
as Only limited variations in domain sizes 

with heating time can therefore be observed in ex¬ 
periments. Amongst inorganic materials, a large 
body of the phase separation literature is devoted 
to silicate melts [Ij. The seminal theoretical work 
of Cahn [19| finds its inspiration in the geometry 
of phase-separated sodium borosilicate glasses. Ki¬ 
netic studies in silicates evidenced diffusive coarsen¬ 
ing both for a droplet |20l[3] and an interconnected 
microstructure |2T|. Because of the limited range 
of sizes that are observable during a evolution, 
a more successful approach for controlling the mi¬ 
crostructure size consists in leveraging thermody¬ 
namics instead of kinetics, by changing the compo¬ 
sition in order to modify the value of the critical 
temperature. Martel et al. jUj, for example, found 
that the typical size of phase-separated droplets de¬ 
creased exponentially with alumina content in cal¬ 
cium aluminosilicate glasses. A similar strategy has 
been applied to metallic glasses |23|. 

However, another coarsening mechanism has 
been proposed for larger domain sizes, namely hy¬ 
drodynamic coarsening, that can occur because of 
viscous flow driven by Laplace pressure [24] . This 
regime is characterized by a linear growth of domain 
sizes with time, that should be much more flexible 
to tune the microstructure, for a given composition. 
Several numerical studies [15| US] have addressed 
hydrodynamic coarsening; however, experimental 
observations of this regime are confined to the 
field of colloid-polymer mixtures [23 [28] and poly¬ 
mers [HES], with no single experimental evidence 
in the field of inorganic materials. Since mecha¬ 
nisms such as crystallization (in alloys and glass¬ 
forming melts), gravity or convection and shear of¬ 
ten disrupt coarsening and result in different mor¬ 
phologies, one may wonder if a linear growth of do¬ 
mains can be observed at all in inorganic materials. 

New insights on phase separation paralleled the 
emergence of experimental techniques in materials 
science. Whereas the first kinetic studies were real¬ 
ized in the Fourier space, using light diffusion tech¬ 
niques EH 11 ], direct access to the morphology of the 
phases is now possible thanks to imaging techniques 
such as electronic microscopy [32] , atomic force mi¬ 
croscopy [21] |20] , X-ray tomography [33] or atom 
probe tomography [M] [23]. In-situ synchrotron 
microtomography [35] is a technique of choice for 


studying microstructure formation E11I37] , since it 
provides both the 3-D microstructure and its topol¬ 
ogy, and successive snapshots of its evolution. Re¬ 
cently EH] ) our team started using in-situ tomogra¬ 
phy to study phase separation in silicate melts. 

In this work, we investigate the coarsening stage 
of a barium borosilicate glass-forming melt, far 
above its glass transition. Thanks to in-situ to¬ 
mography, we evidence a linear growth of the mi¬ 
crostructure with time, characteristic of viscous 
coarsening, that is shown to stem from a series of 
hydrodynamic pinch-offs. In Section]^ we present 
the materials and methods used in this work, in 
particular the characterization of the glasses and 
the in-situ tomography experiments. Section fo¬ 
cuses on the coarsening kinetics during isothermal 
treatments observed by in-situ microtomography, 
for an interconnected geometry. Section studies 
the evolution of the 3-D geometry of the phases 
during the coarsening process. Finally, we discuss 
in Section[^the relative importance of diffusive and 
viscous transport in our system. 

2. Materials and methods 

2.1. Glass composition and preparation 

In this work, we focus on the ternary diagram 
Si 02 — BaO — B 2 O 3 . While alkali borosilicates 
have been much studied in the literature, alkaline- 
earth borosilicates are known to have a higher crit¬ 
ical temperature, and a larger immiscibility do¬ 
main [3] . The immiscibility diagram determined by 
Levin and Cleek El] is shown in Fig.[^ The critical 
temperature of the system is above 1470°C, and a 
large fraction of the immiscibility dome is located 
above liquidus surfaces (not shown here), resulting 
in stable liquid-liquid phase separation. 

The ternary composition studied here is found 
in Table It corresponds to a composition well 
inside the immiscibility domain in order to obtain 
an interconnected microstructure, and to a liquidus 
temperature of 1180°C. 

In order to study submicronic microstructures, 
we have studied a second glass composition with 
a small percentage of alumina [H] [ 22 ] , denoted as 
quaternary glass in Table [l] and subsequently. Since 
the quaternary diagram has not yet been charac¬ 
terized in the literature, we determined the dome 
temperature for this composition by thermal treat¬ 
ments at different temperatures, followed by a rapid 
quench in water. The dome temperature was esti¬ 
mated to be 1175°C (±25°C). 


2 


glass name 

SiOa 

BaO 

B 2 O 3 

AI 2 O 3 

ternary glass 

59 

21.5 

19.5 

0 

Si-rich phase 

80 

2.7 

17.3 

0 

1270-Ba-rich phase 

43 

32 

25 

0 

1070-Ba-rich phase 

32.8 

41.1 

26.1 

0 

quaternary glass 

57 

21 

20 

2 

droplet glass 

66 

14 

20 

0 


Table 1: Composition (in weight percent) of the different barium borosilicate glasses. 


Si02 



Figure 1: Phase diagram of the Si02 — Ba 0 —B 203 system, in 
weight fractions. Isotherms of the immiscibility dome (in or¬ 
ange) are plotted after [39], while the tie-line corresponding 
to composition 1 is drawn from electron microprobe measure¬ 
ments. The composition used in tomography experiments 
corresponds to the ternary glass point. 

For the elaboration of the glasses, batches of 1 kg 
were melted at 1600° C from raw materials (reagent- 
grade barium carbonate and boric acid, and E400 
silica) in a gas-fired lab furnace. Oxygen bub¬ 
bling was used to ensure a good homogeneity of the 
melt. Glasses were quenched in air and annealed at 
650°C. During the quench from 1600° C to room 
temperature, a small microstructure appeared be¬ 
cause of phase separation. The typical size of this 
initial microstructure is 5 /rm for the ternary glass, 
and 100 nm for the quaternary glass (because of a 
much lower dome temperature). 

2.2. Tie-line 

Determining phase compositions in demixed 
glasses is usually a difficult endeavour because of 
the submicronic size of the domains. Typical meth¬ 
ods use selective dissolution in acid m or com¬ 


pare DSC plots to previously-determined glass- 
transition isotherms [H], but such methods often 
lack precision. In our system, large sizes of the 
domains (several tens of microns) can easily be ob¬ 
tained, which made it possible to analyze the com¬ 
positions by electron microprobe. Thermal treat¬ 
ments were realized at 1070° C and 1270° C for the 
ternary glass. The composition of the phases was 
analyzed using a Cameca SX 100 microprobe at 15 
kV and 15 nA. The beam was unfocused with a 
spot size of 20 ^m in order to limit boron migra¬ 
tion. The compositions of the phases are given in 
Table and represented in Fig. Compositions 
are in agreement with the limits of the immiscibil¬ 
ity region proposed by Levin and Cleek [53]. At 
both temperatures, the glass melt separates in a 
silica-rich phase and a barium-rich one, the varia¬ 
tion of the boron content being much smaller. The 
precision of the microprobe measurements - of the 
order of 1 wt% because of the incertitude on boron 
content - is not sufficient to distinguish between 
the silica-rich compositions at 1070°C and 1270°C. 
Interestingly, phases at 1070°C and 1270°C lie on 
the same tie-line, therefore no tie-line rotation is 
observed in such a temperature range (contrary to 
other studies made at lower temperature |43] , closer 
to the glass transition). 

2.3. Physieal properties 

As will be described later, the kinetics of domain 
coarsening is governed by the ratio between the in¬ 
terfacial tension between the two phases, and the 
viscosity of the most viscous phase. In order to 
characterize the physical properties of the separated 
phases, glasses with the composition of the silica- 
rich phase, and the barium-rich phases at 1070°C 
and 1270°C, were melted from raw materials using 
the same procedure as for the ternary glass. 

Viscosity. The viscosities of the barium-rich phases 
were measured using a high-temperature Couette- 
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hydrodynamic simulations. The numerical method 
is described in the Appendix. The order of mag¬ 
nitude is consistent with the scarce literature on 
measurements of interfacial tension between glassy 
melts [33 135] , albeit at the lower end of reported 
values. 

Densities. The density of the separated phases 
were measured from the bulk synthesized glasses, 
using Archimedes principle. With a density differ¬ 
ence of Ap ~ 10^ kg.m“^ between the Si-rich and 
the 1070-Ba-rich phase, we estimate the Bond num¬ 
ber 


Figure 2: (a) Schematic of the falling sphere experiment. 
The diameter of the sphere is 1 cm, and the inner diameter 
of the crucible is 5 cm. (b) Arrhenian plot of measured 
viscosities for the silica-rich phase (♦), 1270°-barium-rich 
phase (■) and 1070°-barium-rich phase (•). 


flow apparatus. For the silica-rich phase, the vis¬ 
cosity was too high to use the Couette set-up. We 
therefore resorted to falling-sphere viscosimetry for 
the silica-rich phase. Fig.[^(a) shows a schematic of 
the experiment. A transparent silica crucible con¬ 
taining the glass melt is placed inside an electric 
furnace with a transparent silica window. A plat¬ 
inum sphere is tracked with a camera while it falls 
through the melt, and Stokes’ law is used in order 
to retrieve the viscosity of the melt. Corrections 
were applied to account for the finite size of the 
crucible [ 33 ] and the accuracy of the method was 
checked using a soda-lime glass of known viscosity. 

The viscosities of the different phases are repre¬ 
sented on an Arrhenian plot in Fig. (b). A very 
important viscosity contrast is observed between 
the silica-rich phase and the barium-rich phase, up 
to five orders of magnitude. 

In order to supplement viscosity data, dilato- 
metric measurement were realized on the silica-rich 
phase and the 1070°C-barium-rich phase. The glass 
transition and the dilatometric softening point (cor¬ 
responding to a viscosity of about lO'^'^Pa.s) were 
obtained with this method, with values of respec¬ 
tively 570 and 690°C for the silica-rich phase, and 
640 and 670°C for the 1070°-barium-rich phase. 
Such measures confirm the large extent of the tem¬ 
perature range at which phase separation occurs. 

Interfacial tension. We estimated the value of the 
interfacial tension 7 between the two phases to be 
of the order of 10“^ J.m"^, from the comparison of 
experimental images of the phases and numerical 


ApgP 2 

Bo= -AAl— ~ 10-2 


( 1 ) 


for a typical size £ of 100 ^m, which is the maxi¬ 
mum size that is reached in our coarsening exper¬ 
iments. The Bond number measures the ratio be¬ 
tween gravitational and interfacial forces. Its small 
value means that gravity can be neglected in com¬ 
parison to Laplace pressure. 


2 . 4 . Imaging methods and thermal treatments 

In-situ X-ray microtomography experiments were 
realized at the ID19 beamline of the ESRF syn¬ 
chrotron. The imaged samples were 2-mm-diameter 
cylinders machined from the bulk ternary glass and 
put inside 3-mm-outer-diameter alumina crucibles. 
On the beamline, we used the dedicated “Ecole 
des Mines” furnace, that can heat millimetric-sized 
samples to temperatures up to 1500°C. During in- 
situ experiments. X-rays can get through the fur¬ 
nace thanks to a thin silica window. Samples are 
glued on top of an alumina rod fixed on a Leu¬ 
ven fast rotation stage, and are allowed to rotate 
freely inside the fixed furnace thanks to a hole at 
the bottom of the furnace. At the beginning of a 
thermal treatment, the sample is introduced into 
the already hot furnace, ensuring a fast heating 
rate: we measured that it only takes a few tens 
of seconds for the sample to reach working temper¬ 
atures above 1000°C. Isothermal experiments were 
realized at temperatures ranging from 1030° C to 
1330°C. 

Pink beam with an energy of 32 keV is used to 
achieve sufficient transmission. At this energy, the 
absorption coefficients of the two phases (resp. 0.08 
mm“^ and 1 mm“^ for the silica-rich and the 1070°- 
barium-rich phases) are very different, resulting in 
a good absorption contrast that facilitates further 
image processing. The X-ray photons transmitted 
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through the sample were converted to visible light 
by a 25 /rm-thick LuAG scintillator. The pixel size 
is 1.1/im. The PCO Dimax camera is equipped 
with a fast CMOS detector; acquisition times for 
one radiography varied from 6 to 30 ms depending 
on the working temperature, resulting in total ac¬ 
quisition times from 4 to 10 s for the total number 
of radiographies acquired over a rotation of 180°. 

3-D absorption images were reconstructed from 
the set of radiographies using a standard fil¬ 
tered back-projection algorithm [17]. Python’s 
scikit-image [48] was used for further processing 
of 3-D images. Voxels were attributed to the two 
phases of the glass using the following segmentation 
procedure: a first denoising step is realized with a 
total-variation filter |1S|; then the Random Walker 
algorithm [^ is used to segment the two phases. 
3-D visualizations of the minority-phase surface are 
realized with the Mayavi [51] Python package. 

For complementary experiments at scales out of 
the reach of parallel-beam synchrotron tomography, 
thermal treatments were realized on centimetric- 
sized samples in an electric furnace. The samples 
were quenched and impregnated in an epoxy resin. 
Cross-sections were polished and imaged with a 
Zeiss DSM 982 Gemini scanning electron micro¬ 
scope, operated at 20 kV. 

3. Viscous coarsening 

The typical evolution of a sample is shown on 
Fig-i for a temperature of 1230°C. The surface 
of the barium-rich phase is represented on Fig. 
with colors encoding the value of the local mean 
curvature (from dark green to bright yellow for in¬ 
creasing curvature). The same cubic sub-volume of 
linear size 250 /im is represented at different times 
during an in-situ experiment, from 80 to 480 s. In 
this experiment, the volume fraction of the barium- 
rich phase is measured to be 49 ± 1%, so that it is 
in slight minority. 

An interconnected microstructure is observed for 
both phases. The minority barium-rich phase is 
composed of smooth threads of fluid connected to¬ 
gether by thicker nodes. Such an interconnected 
microstructure is often attributed to spinodal de¬ 
composition, however our observations only start 
when the compositions of the phases are fixed. A 
few disconnected barium-rich domains are observed 
as well (and represented in purple in Fig. [^, but 
they represent a very small volume fraction of the 



Figure 3: Tomography images of domain coarsening, for an 
experiment at 1230°C and heating times of 80, 160, 320 and 
480 s. The size of the square box is 250 /im. Only the 
minority (barium-rich) phase is represented. Shades of green 
encode the local mean curvature, from dark green to yellow. 
Isolated domains are identified in purple. 

phase (always less than 1 % throughout the experi¬ 
ment). No disconnected domains are ever observed 
for the silica-rich phase. As time increases, we ob¬ 
serve the coarsening of the typical size of the mi¬ 
crostructure (Fig. 1^, while the morphology of the 
phases continues to look quite similar to the initial 
state, up to a rescaling of the typical size. 

We measure the typical size £ of the barium-rich 
phase as its volume V to surface S ratio: 


so that i is equal to the radius for a sphere. The 
surface area is measured from the binarized image 
thanks to local histograms of voxel configurations, 
as proposed by Lang et al. [^. The segmentation 
process often results in surfaces with a small-scale 
roughness, at a scale below the true spatial reso¬ 
lution of the experiment. Therefore, the surface 
is smoothened with a Gaussian filter of standard 
deviation of 1 pixel before the surface area is com¬ 
puted. This smoothing step reduces the sensitivity 
of i to the measure and segmentation noise. We 
have compared £ to other measurements of the aver¬ 
age local scale derived from the correlation function 
of the binarized image, or the chord length distri¬ 
bution IHIMKIl- All quantities resulted in similar 
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Figure 4: Evolution of the characteristic length £(t). Main 
axes: growth of i measured in in-situ tomography experi¬ 
ments at several temperatures. Inset: growth of £^^^1 mea¬ 
sured post-mortem thanks to SEM imaging. 

results [55] (up to constant prefactors of order one, 
coming from the different definitions). 

The evolution of £ with time is shown on Fig. [^ 
for experiments realized at different temperatures, 
from 1180°C to 1280°C. A linear evolution with 
time is observed, up to typical sizes as high as 
80 fiia. Since the initial microstructure has a size of 
about 5 /im, the typical size is multiplied by more 
than one order of magnitude during the duration of 
an experiment. The coarsening rate increases with 
temperature, and several experiments performed at 
the same temperature show that the value of the 
coarsening rate is well reproduced (see Fig. |^. In 
order to characterize the temperature-dependence 
of the coarsening rate, we have plotted £ in an Ar- 
rhenian diagram in Fig. for 11 different experi¬ 
ments realized from 1030 to 1330°C. We observe an 
Arrhenian temperature dependence, with an acti¬ 
vation energy of the order of 300kJ.mol“^. Since 
the immiscibility dome widens towards lower tem¬ 
peratures, experiments at different temperatures in 
Fig. [5] correspond to different volume fractions of 
the barium-rich phase, from 39 to 52 %. Such vari¬ 
ability does not seem to affect the Arrhenian be¬ 
havior. 

In order to test the validity of the linear coarsen¬ 
ing regime at smaller scales, we measured the typi¬ 
cal size of interconnected phases for thermal treat¬ 
ments at 1000° C of the quaternary glass with alu¬ 
mina, using post-mortem SEM imaging. Since £ 
cannot be directly measured from 2D images, we 
measured as a proxy for £ the average chord size 
the chord size distribution being a stereo- 



Figure 5: Compararison between experimental and theoret¬ 
ical coarsening rate, plotted in an Arrhenian diagram. Cir¬ 
cles: linear growth rates v = £ measured in experiments at 
temperatures from 1030 to 1330°C. Diamonds: ratio ■y/rj, 
determined using the value of 7 obtained numerically (see 
Appendix) and measured viscosities. 


logical measure that should be identical in 2D and 
3D for statistically isotropic systems. As shown in 
the inset of Fig. a linear growth of £^Yy with time 
is observed for this set of experiments as well, for 
values of £^^1 ranging from 1 to 3 fim. The first data 
point, corresponding to a scale of 300 nm, does not 
match as well the linear fit. This slight discrepancy 
at short times and small scales is likely to come 
from variations in the initial size of the microstruc¬ 
ture, since our post-mortem observations were real¬ 
ized on different samples. We conclude that a linear 
coarsening can be observed at scales ranging from 
1 to 80 /im, including scales below the resolution of 
our tomography experiments. 

The observed linear evolution of £ suggests that 
viscous coarsening dominates over diffusive mecha¬ 
nisms. This may come as a surprise, since diffusive 
coarsening has been the only mechanism reported 
so far in silicate melts jl). 3-D visualizations of 
the evolution of interfaces (Fig.[^ unravel the local 
mechanisms resulting in coarsening. Inside curved 
liquid bridges, Laplace pressure causes fluid to flow 
from regions with a high mean curvature towards 
regions with smaller curvatures. This results in a 
thinning of some liquid bridges (Fig. § (b)), up to 
the break-up of the liquid filament (Fig. §(c)). Af¬ 
ter the break-up, fluid in the broken ends retracts 
because of capillarity and therefore flows in adja- 
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Figure 6: Close-up visualization of a breaking liquid bridge. 
Surfaces of the percolating domain are color-coded according 
to mean curvature, from green to yellow. A few disconnected 
droplets are observed as well (in purple). 3-D images of 
linear size 140 microns, taken at times 408, 432, 440, and 
464 s during an experiment at 1280°C. 


cent regions (Fig.[^(c) and (d)), resulting in a local 
coarsening. Local break-up events of viscous liquid 
bridges are also known as pinch-off in the fluid me¬ 
chanics literature |56] . They were observed using 
confocal microscopy in a colloid-polymer mixture 
by Aarts and collaborators m- Further evidence 
for the weak influence of diffusion on coarsening is 
given by the fact that a few isolated domains (dis¬ 
played in purple in Fig. keep a constant volume 
through time, so that evaporation can be considered 
negligible up to the measurement precision. Visual 
observation also sheds light on the origin of the few 
isolated domains, that break up from the contin¬ 
uous phase during the retraction of long terminal 
filaments [55] (see also Fig. [^. 

Thanks to such direct observations of coarsening 
mechanisms, it is possible to revisit the theoretical 
argument of Siggia [24] for deriving a linear coars¬ 
ening rate with time. One balances the dominant 
terms of the Navier-Stokes equation, that are the 
Laplace pressure gradient and the viscous dissipa¬ 
tion, when inertia is negligible: 

Pbaplace T 

with 77 the largest viscosity of the two phases, 7 the 
interfacial tension, v a typical velocity of the fluid 


and the Laplace pressure pLapiace ~ 7 /L Siggia ar¬ 
gues that V ~ d£/dt, since £ is the typical scale of 
fluid domains. In Fig. [^ v would be a typical ve¬ 
locity at which fluid flows from the bridge to the 
node, either when the bridge is still connected or 
when the broken ends retract. The growth of £ in 
neighboring liquid bridges is due to the redistribu¬ 
tion of the fluid from the thinning, then breaking, 
liquid bridge. Replacing v, one obtains the equation 
proposed by Siggia: 

£{t) ^ U (4) 

In order to test whether Eq. describes well the 
coarsening dynamics observed here, we have com¬ 
pared in Fig. 1^ the evolution of the linear prefactor 
with the evolution with temperature of the ratio 
7 / 77 . In Fig. the linear prefactors are plotted 
in an Arrhenian diagram, for temperatures ranging 
from 1030 to 1330°C. As for the viscosity, the kinet¬ 
ics of fluid rearrangement can be considered to be 
governed by the largest viscosity. The estimated in¬ 
terfacial tension over the measured viscosity of the 
silica-rich phase is therefore superimposed on Fig.[^ 
We observe a very good agreement between £{t) and 
7 / 77 : the activation energies of the two quantities 
are very close. Moreover, the proportionality coef¬ 
ficient between £{t) and 7/77 seems to be close to 
one. The good quantitative agreement between our 
measurements and Eq. 0 is a convincing evidence 
that viscous coarsening is a mechanism present in 
phase-separated silicate melts. 


4. Self-similarity of the microstructure 

An important hypothesis underlying theoretical 
descriptions of coarsening is that the geometry of 
the phases remains statistically the same up to a 
rescaling of the typical size £ - a hypothesis called 
dynamic scaling m- In particular, this argument 
is used by Siggia |53| to justify that £ is equal to 
the velocity induced by the Laplace pressure. 

Statistical characterizations of the domain geom¬ 
etry allow us to test the dynamic scaling hypoth¬ 
esis. In order to characterize the geometry of the 
system, we compute histograms of surface curva¬ 
tures, for an experiment at 1180°C. Distributions 
of curvatures are a common tool for a local char¬ 
acterization of interfaces, in particular for coarsen¬ 
ing studies im |59|. Principal curvatures ki and 
K 2 are computed by fitting a quadric surface to the 
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Figure 7: Histograms of mean and Gaussian curvatures, for 
an experiment at 1180°C, observed at thermal treatment 
times of 15, 28, 48 and 60 min, from dark blue to light blue, 
(a) Mean curvature H. (b) Renormalized mean curvature 
Hi. (c) Gaussian curvature K. (d) Renormalized Gaussian 
curvature Ki^. 


neighbourhood of points lying on the interfacial sur- 
face [60l [61]. Mean and Gaussian curvatures are 
defined respectively by 

H= -^-; K = kiK2. (5) 

Histograms of H and K measured at times from 15 
to 60 minutes are plotted in Fig. Histograms of 
curvatures get narrower with time because of global 
coarsening, as the characteristic length £ grows from 
25 to 75 microns from the first to the last image 
used to build the histograms. The range of curva¬ 
ture values reflects the different possible morpholo¬ 
gies of the interfacial surface, with rounded tips 
{K > 0, H < 0), tubes {\K\ < ^"2, H < 0) or 
saddles {K < 0). For both H and K, the peak 
of histograms correspond to tubular morphologies, 
and the tails to either tips of such tubes, or junc¬ 
tions between them. Fig. (b) and (d) shows the 
histograms of mean and Gaussian curvatures renor¬ 
malized by the characteristic length £, Hi and Ki"^. 
A good collapse on a master curve is observed both 
for mean and Gaussian curvature, demonstrating 
that, statistically speaking, the morphology of the 
phases does not change over time, except for their 
average size. The validity of dynamic scaling means 
that it is possible to elaborate morphologies that 
are statistically similar, but with sizes that are dif¬ 
ferent. Such property can be interesting for appli¬ 
cations but also for designing model materials. 

Nevertheless, one should mention that deviations 
from dynamic scaling can be observed at longer 
times, because of the progressive fragmentation 
of disconnected droplets from the percolating low- 
viscosity phase. Such domains can be observed in 
Fig. [3] and We have studied deviations from 
dynamic scaling in a previous paper [35]. Recent 
work [62] has shown that the fragmentation rate in¬ 
creases drastically when the volume fraction of the 
low-density phase decreases. In the set of experi¬ 
ments presented here, the volume fraction is large 
enough that a small fragmentation rate hardly dis¬ 
rupts dynamic scaling for the duration of experi¬ 
ments. 


5. Discussion 

Our tomography experiments demonstrate that 
viscous flow induced by Laplace pressure is the pre¬ 
vailing coarsening mechanism over a wide range of 
lengthscales. Then, why did not previous studies 
on phase separation in glass-forming melts [M1I2I] 
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observe a linear-with-time coarsening? In order to 
investigate the relative effects of diffusion-induced 
and viscous-ffow-induced coarsening, we carried on 
a complementary set of experiments where diffu¬ 
sion was the sole coarsening mechanism. To do so, 
we synthesized a composition on the same tie-line 
as the ternary glass, but sufficiently off-centered 
towards the silica-rich phase so that the compo¬ 
sition lies inside the binodal domain (see Fig. [^, 
and is characterized by a droplet morphology (see 
Fig. § (a)). After the quench from the stable re¬ 
gion to below the glass transition, a microstructure 
of barium-rich droplets inside a silica-rich matrix 
is visible. In the absence of an interconnected mi¬ 
nority phase, the prevailing transport mechanism 
is chemical diffusion through the boundary of the 
domains (droplets coalescence because of Brownian 
motion [64] is bound to be slower because of the 
micronic size of droplets and the very large viscos¬ 
ity of the matrix). Thermal treatments were real¬ 
ized for cm-sized glass cylinders at 1250° C for times 
ranging from 20 minutes to 6 hours. Polished sur¬ 
faces cut through the samples were observed post¬ 
mortem using SEM imaging. Typical SEM images 
are shown in Eig. (a). The histogram of areas of 
the barium-rich droplets was converted to an his¬ 
togram of volumes of 3-D spheres, using a classical 
stereologic computation j^. The evolution of the 
average domain size with time is shown in Fig.[^(b); 
it is consistent with a evolution. In addition to 
the different kinetic law, we note that the typical 
size of droplets is much smaller than for the inter¬ 
connected structure at a similar temperature (see 
Eig. 1^. This size discrepancy confirms that, when 
topology allows it, viscous flow is a much more ef¬ 
ficient mechanism than diffusion for our system. 

For a diffusion-controlled coarsening, diffusion 
models predict that the growth of i evolves as m 


( 0 ) 

where D is the effective molecular diffusivity, D 
a molecular volume, k the Boltzmann constant 
and T the temperature. Fitting the size values 
of Fig. (b) with Eq. we obtain an esti¬ 

mate for the product DD. Choosing the molecu¬ 
lar volume D ~ gives D ~ 2.10“^^ m^s“^. 

The theoretical crossover length £* between diffu¬ 
sive growth and hydrodynamic growth is commonly 
estimated by equating the interfacial velocity in the 
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Figure 8: (a) SEM images of coarsening for a droplet mor¬ 
phology, for thermal treatments of 0, 20, 120 and 360 min at 
1250°C (left to right), (b) Evolution of the average droplet 
radius with 


two regimes [66], resulting in 



(7) 


Eor our system at 1250°C, the corresponding nu¬ 
merical values is I* ~ 10 nm, much below the ini¬ 
tial microstructure size in tomography experiments. 
The corresponding crossover time t* = 
would be of the order of 0.1 s. Observing only the 
hydrodynamic regime at micronic scales was there¬ 
fore to be expected. 

In this respect, it is instructive to compare 
such orders of magnitude to the ones of another 
study [53] where both diffusivity and viscosity were 
estimated. Simmons et al. [53] measured the evo¬ 
lution of the macroscopic viscosity resulting from 
coarsening in sodium borosilicate glasses. Using the 
values obtained for rj and D in [53] and Eq. 0 , one 
obtains a crossover scale of a few microns and an 
hypothetical crossover time of thousands of days. 
Observing the hydrodynamic regime requires t* to 
be accessible in an experiment, and hence the vis¬ 
cosity of the glass not to be too large. Since most 
previous studies in glass-forming melts were real¬ 
ized slightly above the glass transition, the viscosity 
was always very high (> lO^^ Pa.s) and the hydro- 
dynamic regime was out of reach. An interesting 
perspective of our study would be the investigation 
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of the crossover of the diffusive and hydrodynamic 
regimes. Such investigations would require to con¬ 
trol the critical temperature so that the viscosity 
just below Tc would result in a crossover time well 
resolved experimentally. 

Let us finally note that simple models relat¬ 
ing the viscosity and the diffusivity, such as the 
Stokes-Einstein relation, or the Eyring model m 
for silicate glasses, result in a crossover length 
I* of the same order as a typical molecular size. 
Observations of diffusive coarsening over larger 
scales [211 [63] therefore hint at more complex dif¬ 
fusive mechanisms that are decoupled from viscos¬ 
ity - a decoupling that is well known for silicate 
glasses close to their glass transition [^ [51] and 
that has been recently discussed for metallic glasses 
as well [70]. 

6. Conclusions 

Using in-situ microtomography, we have shown 
that viscous flow prevails as the dominant coars¬ 
ening mechanism in phase-separated silicate melts, 
for temperatures far above the glass transition. A 
linear growth of the typical size has been observed 
for a large range of scales, from 1 to 80 p,m. We 
have verified that the morphology remains statisti¬ 
cally self-similar during coarsening, up to a rescal¬ 
ing by the characteristic size. Whether these mech¬ 
anisms apply to other inorganic materials, such as 
phase-separated bulk metallic glasses [3 [71] , is an 
open question of interest for the fundamental un¬ 
derstanding of such materials. 

The transition between the diffusive and vis¬ 
cous regime is an interesting question open by this 
work, that will require to understand which dif¬ 
fusive mechanisms control the coarsening at small 
scales. Future work will focus on smaller volume 
fractions on the less-viscous phase, for which more 
fragmentation of droplets is observed. 
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Appendix: estimation of interfacial tension 

In order to estimate the interfacial tension be¬ 
tween the two phases, we performed numerical sim¬ 
ulations with the Gerris software [73173], a finite- 
volume solver of Navier-Stokes equations. As an 
initial condition, a mesh was built from the seg¬ 
mentation of a single domain in an experiment at 
1180°C (see Fig. |^. Physical parameters of the 
simulation were set to one, including the size of the 
box, the interfacial tension between the phases and 
the viscosity of the most viscous phase. The vis¬ 
cosity of the least viscous phase was set to 10 “^, 
ensuring a high viscosity ratio between the phases. 
The viscosity ratio in the simulation is not as large 
as the experimental one, since too large a dis¬ 
crepancy between the two viscosities leads to pro¬ 
hibitive computational times. Such parameters en¬ 
sure that the Reynolds number is smaller than 0.1. 
For Reynolds numbers much smaller than 1, in the 
Stokes flow regime, inertia is negligible and time 
is only a parameter, with boundary conditions de¬ 
termining completely the further evolution of the 
system. 

Fig. [9] shows the evolution of the initial condi¬ 
tion in the experiment, as well as in the simulation 
of Navier-Stokes equations. A striking similarity 
is observed between the experimental and numer¬ 
ical geometries. In particular, the domain breaks 
up at the same place, while the lower end of the 
domain retracts in both cases. Such observations 
confirm that hydrodynamic flow is the mechanism 
responsible for the observed coarsening. Moreover, 
the proportionality factor needed to adjust experi¬ 
mental and non-dimensional numerical times gives 
the value of the interfacial tension in the experi¬ 
ment, that is found to be 7 = 10“^J.m“^. We 
suppose that the interfacial tension does not vary 
much with temperature, which is a reasonable as¬ 
sumption compared to the much greater variations 
of the viscosities with temperature. 
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